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Abstract 

An improved composite-boson theory of quantum Hall ferromagnets is pro- 
posed. It is tightly related with the microscopic wave-function theory. The 
characteristic feature is that the field operator describes solely the physical de- 
grees of freedom representing the deviation from the ground state. It presents 
a powerful tool to analyze excited states within the lowest Landau level. Exci- 
tations include a Goldstone mode and nonlocal topological solitons. Solitons 
are vortices and Skyrmions carrying the U(l) and SU(2) topological charges, 
respectively. Their classical configurations are derived from their microscopic 
wave functions. The activation energy of one Skyrmion is estimated, which 
explains experimental data remarkably well. 

I. INTRODUCTION 

The quantum Hall (QH) effect IJ is a remarkable macroscopic quantum phenomenon ob- 
served in the two-dimensional electron system at low temperature T and in strong magnetic 
field B. The Hall conductivity is quantized with extreme accuracy and develops a series of 
plateaux at magic values of the filling factor v = 27ihp /eB. Here, p is the average electron 
density. The QH effect comes from the realization of an incompressible quantum fluid. The 
composite-boson (CB) picture p|d4| and the composite-fermion picture |||6| have proved to 
be quite useful to understand all essential aspects of QH effects. Electrons may condense 
into an incompressible quantum Hall liquid as composite bosons. The QH state is such a 
condensate of composite bosons, where quasiparticles are vortices |7[]. 

When the spin degree of freedom is incorporated, a quantum coherence develops spon- 
taneously provided the Zeeman effect is sufficiently small, turning the QH system into a 
QH ferromagnet. New excitations are Skyrmions. Skyrmions were initially considered as 
solutions of the effective nonlinear sigma model ||, and later studied also in a Hartree-Fock 
approximation M. Their existence has been confirmed experimentally [|T0| , |TT| 1 . 



In this paper we present a field-theoretical formulation of QH ferromagnets and 
Skyrmions based on an improved composite-boson (CB) theory [12|], proposed based on 
a suggestion due to Girvin QTJ] , Read || and Rajaraman et al. [|J. In this scheme, the 
field operator describes solely the physical degrees of freedom representing the deviation 
from the ground state, and the semiclassical property of excitations is determined directly 
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by their microscopic wave functions. This clarifies the relation between the microscopic 
wave function and the classical Skyrmion, which is missing in the effective theory ||. It 
is explicitly shown that the Skyrmion excitation is reduced to the vortex excitation in the 
small-size limit. We estimate its activation energy, which accounts for the observed data 
due to Schmeller et al. |IT| remarkably well provided certain physical assumptions are made. 
Thus, our field-theoretical approach presents a new insight into the Skyrmion physics. 

This paper is composed as follows. In Section |TJ we give an overview of the improved CB 
theory. In Section |III|CB fields are defined. First, as in the standard CB theory JJ, we attach 
an odd number of flux quanta to an electron by way of a singular phase transformation. We 
call the resulting electron-flux composite the bare composite boson. In order to soften the 
singularity brought in, we dress it with a cloud of an effective magnetic field that bare 
composite bosons feel. The resulting object turns out to be the dressed composite boson. 
In Section |V| the relation is established between the electron wave function and the CB 
wave function. We also verify that the ground state is given by the Laughlin state within 
the semiclassical approximation. In Section [V| we make a semiclassical analysis of vortex 
excitations. In Section [VT] we include the spin degree of freedom and show that a quantum 
coherence develops spontaneously when the Zeeman effect is sufficiently small. In Section |VII 
we discuss vortex and Skyrmion excitations in the QH ferro magnet. The Skyrmion classical 
configuration is also derived directly from its microscopic wave function. We evaluate the 
excitation energy of one Skyrmion and compare it with experimental data. Throughout the 
paper we use the natural units K = c = 1. 



II. OVERVIEW 

We start with a review of the improved CB theory applyied to the spin-frozen QH system 
12 1 . We denote the electron field by ip(x) and its position by the complex coordinate 
normalized as z = (x + iy) /2£b- Any state |$) in the lowest Landau level at the filling factor 
v — 1/m (m odd) is represented by the wave function, 

$[x] = {0\if>(x 1 )---il>(x N )\$)=u[z]$ Ijli [x], (2.1) 

where 

N 

$lnM = - ^) m exp[- ^ \z r \ 2 ] (2.2) 

r<s r=l 

is the Laughlin function, and u[z] = ou(zi, Z2, ■ ■ ■ , zjv) is an analytic function symmetric in 
all N variables. The mapping from the fermionic wave function to the bosonic function 
u>[z] defines a bosonization. We call the underlying boson the dressed composite boson and 
denote its field operator by <p(x). The field operator turns out to be the one considered first 
by Read H and revived recently by Rajaraman et al. [[|. We derive that 

%[x] = (0\<p{x 1 )~-<p(x N )\$)=w[z]. (2.3) 

The Laughlin state is represented by $tp[x] = 1. 
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A. Vortices 

A typical vortex state is represented by &<p[x\ = Yl r z r- When the wave function is 
factorizable, io[z\ = Ylr u { z r)i as i n this example, we obtain a semiclassical equation, 

(<p{x))=u(z). (2.4) 

This is a highly nontrivial constraint since it connects directly the classical field (<p(x)) to 
the microscopic wave function u(z). It dictates how the electron density p(x) is modulated 
around the topological excitation. 

The topological charge density Q(x) is determined in terms of the classical field. It is 
the vorticity, Q(x) = Q v (x) given by 

Q V (x) = ^ijdidj H<pW) = (2.5) 

for a vortex sitting at the origin. The topological soliton induces the density modulation 
according to the soliton equation, 

^V 2 \np(x)-p(x)+ Po = uQ(x), (2.6) 

as follows from the semiclassical constraint ( |2.4| ). The topological soliton carries the electric 
charge uqe, where q = J d 2 xQ(x) is the topological charge. 



B. Skyrmions 

We next summarize the idea of an improved CB theory applied to QH ferromagnets. We 
denote the spin component by the index a(=f,|). Any state at v — 1/m is represented 
by the wave function similar to fl2.1|) . Let us explicitly consider the case when the spinor 
component is factorizable, 

* ,pta M=n(:![^) r ^N[4 < 2 - 7 > 

The ground state is given by 

$f n [x]=n(j) <W4 (2.8) 

From a set of two analytic functions u a (z) we construct the complex-projective field n{x) 
whose two components n a (x) are 

n a (x) = 1^2^ (2 9) 

This represents the most general Skyrmion configuration [13 together with the asymptotic 
boundary condition = 1 and v} = 0. The normalized spin field, or the nonlinear sigma 
field, is defined by 
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s a (x) = n(x)^T a n(x), 



(2.10) 



with r a the Pauli matrices. The Skyrmion configuration is indexed by the Pontryagin number 
13] , whose density is 



Q p (x) = ^-e abc e ij s a (x)d i s b (x)d^ c (x). (2.11) 
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The Skyrmion excitation modulates not only the spin density but also the electron density. 
The electron density modulation is determined by the same soliton equation as (|2.6|) , where 
Q(x) = Q p {x) is the Pontryagin number density ( |2.11| ). 

The simplest excitation is given by one Skyrmion with scale k sitting at x — 0, whose 
wave function is [ITJJ] 

$ sCW = n( K %) $ ln[x]. (2.12) 

The scale k is to be fixed dynamically to minimize the excitation energy. The Pontryagin 
number density ( [2. Ill ) is calculated for the simplest Skyrmion ( |2.12|) as 



W = Iwr&f- (2 - 13) 

We study the limit k — > 0, where the Skyrmion wave function ( |2.12|) is reduced to the 
vortex wave function in the spin-polarized ground state, and the Pontryagin number density 
is reduced to the vorticity ( |2.5| ), Q p (x) — > S(x), as n — > 0. Hence, the soliton equation 
( |2.6|) yields a vortex configuration for the zero-size Skyrmion. This makes clear the relation 
between the Skyrmion and the vortex in QH ferromagnets. 



III. BOSONIZATION 

The microscopic Hamiltonian for spinless planar electrons in external magnetic field 
(0, 0, —B) is given by 

H = ^mJ dVWft 2 + Py)Hx) + H c (3.1a) 

= M J d2 ^ t ( x )( p - - ip y^ p * + ip y)^( x ) + y Wc + Hc > (3 ' lb) 

where ip(x) is the electron field; Pj = —idj + ev4° xt is the covariant momentum with 
Aj Xt = ^EjkXkB; N is the electron number and uj c is the cyclotron frequency. The Coulomb 
interaction term is 

Ha = ff ^y^M, (3.2) 
2e J \x-y\ 

where g(x) = p(x) — po stands for the deviation of the electron density p(x) = ip^ {x)ip{x) 
from its average value po- 
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The state |<3>) in the lowest Landau level obeys 



(P x + iP y )^(x)m = -±(z + —)l,{x)\Q) = a, (3.3) 

upon which the kinetic Hamiltonian is trivial. We call it the LLL condition. The generic 
solution of this equation yields the iV-body wave function ( |2.1| ) for the state |<3>). We are 
concerned about the state |$) in the lowest Landau level. 

The bosonization scheme was pioneered by Girvin and MacDonald 0. We define the 
field 4>(x) by an operator phase transformation, 

(j)( x ) = e - i0(x V(^), (3.4) 

where Q(x) is the phase field, 

e W .n./dV(.-^). (3.5) 



with the azimuthal angle 6{x — y). When m is an odd integer, we can prove that <f>(x) is 
a bosonic operator. Let us call the underlying boson the bare composite boson. It is a 
hardcore boson satisfying the exclusion principle, 0(x) 2 = 0. The LLL condition (|3.3| ) reads 

(P x + iP y )(f>(x)\$) = 0, (3.6) 

where Pj = Pj + djQ(x) is the covariant momentum for the bare CB field. By solving 
this condition, the iV-body wave function is found to be $$[x] = u;[z]|$lnNI- Namely, the 
standard bosonization is a mapping from the fermionic wave function <f>[x] to a bosonic wave 
function ^[sc], 

$[x] »$4,[x] =co[z]\^[x)\. (3.7) 

The mapping is to attach m Dirac flux quanta mipr, to each electron by way of the phase 
transformation (|3.4|) , where 4>d = 2n/e. Thus, the ground state of bare composite bosons is 
described by the modulus of the Laughlin function, |<&ln[x]|. 

We wish to introduce another CB field (p(x) which makes the LLL condition as simple 
as possible. We set 

<p(x) = e-^'-^VW, (3.8) 

together with an operator A(x) to be determined later. By substituting this into the LLL 
condition (|3.3|) , provided A(x) satisfies 

djA{x) = s jk {d k e(x) + eA c k xt (x)} = e jk d k e{x) - ^ Xj , (3.9) 
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it is reduced to a simple formula, 



(r x + iv y Mx)m = -i-JLipixm = o, (3.io) 
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where Vj is the covariant momentum for the new CB field <p(x): See ( |4.3|) . We solve ( |3.9| ) 

as 



A(x) = mj d 2 y\n{^^-)p(y) - \z\ 2 . (3.11) 

The iV-body wave function ^[x] of dressed composite bosons is obtained by solving the 
LLL condition (|3.10| ), and it is given simply by an analytic function u[z] as in ( [2.3|) . We 
shall verify in the next section that one analytic function u[z] characterizes one state |$) as 
in ( |2.1|) in terms of electrons, or as in ( |3.7| ) in terms of bare composite bosons, or as in ( |2.3| ) 
in terms of dressed composite bosons. 

The effective magnetic potential for the bare composite boson is A| xt + (l/e)(9fcO, which is 
rewritten as (l/e)e k jdjA(x). Bare composite bosons feel the effective magnetic field B e s(x), 

B cS (x) = e~ l V 2 A(x) = m<p D p(x) - B. (3.12) 

The effective field vanishes, (B e g) = 0, on the homogeneous state (p(x)) = po if mcj)DPo = B. 
It occurs at the filling factor v = p^^jB = 1/m. At v — 1/m we rewrite the effective 
magnetic field as 

eB eS {x) = V 2 A{x) = 2irmQ{x), (3.13) 
with g(y) = p(y) — p . It is solved as 

A(x) =mj d 2 y\n(^^-)g(y). (3.14) 

This formula is equivalent to ( p.ll| ) at v — 1/m up to an irrelevant integration constant. 
Since the field <p(x) is constructed by dressing a cloud of the effective magnetic field, we 
have termed it the dressed CB field. 

IV. QUANTUM HALL STATES 

We derive the Hamiltonian in terms of dressed composite bosons by substituting 
into ( |3.1b| ), 



H = 1 M I d2 V(z)(P* - i^vXP* + WvMx) + juj c + H c , (4.1) 
where we have defined 

V t (x) = ^(x)e 2A{x \ (4.2) 

with which p(x) = ip j; (x)il!(x) = (p^(x)(p(x). The covariant momentum Vj = Pj + djQ — idjA 
reads, 

Vj = -idj + e(5 jk - ie jk )Ak, A k (x) = --e kj djA(x). (4.3) 
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As we noticed in ( |3.10|) , it yields 

i d 



with which the Hamiltonian (14.11) is rewritten as 



It is manifestly hermitian. 
The Lagrangian density is 

L = ^{id t - eAf*)ij) -H = <p*(idt - e4f - d t & + id t A)<p - H, (4.6) 

where 7i is the Hamiltonian density and = (Af xt , A^) is the potential of the external 
electromagnetic field. Because the canonical conjugate of <p(x) is not i<p^(x) but i(p^(x), the 
equal-time canonical commutation relations are 

[<p(x), ^{y)\ = 5(x - y), [<p(x), <p(y)} = [^(x), </?%)] = 0. (4.7) 

They are also derived Q by an explicit calculation from those of the electron fields i/j(x) 
and i]r(x) based on the definition ( |3.8| ). 
The CB wave function is defined by 

%[x] = {0\<p(x 1 )<p{x 2 )---<p{x N )\$). (4.8) 



The LLL condition (p. 10 ) implies that the wave function <& v [x] is an analytic function, 



*&ip[x] = u[z]. With the use of the formula ( 3.11]) it is an easy exercise to derive the following 
relation ||, 

V^M^M ■ ■ -<pHxn)\0) = $LNNV t (x 1 )V t (^2) • • ^(x N )\0) t (4.9) 

where $lnN is the Laughlin function ( |2.2|) . 

Because of the commutation relations ( |4.7|) the state |$) associated with the CB wave 
function (|4.8j) is given by 



|$)= l[dx]%[x]^(x 1 )^(x 2 )--- ( pi(x N )\0) (4.10) 
[dx}u[z]^[x}^(x 1 )^(x 2 ) ■ ■■^(x N )\0), (4.11) 



where use was made of Q4.9|), and [dx] = d 2 a:id 2 a;2 • • -d 2 XN- It follows from ( f4.11p that the 



electron wave function is <&[x] = oj[z\$l,n[x], as verifies the basic formulas (|2.1|) 

The semiclassical ground state is the one that minimizes the total energy (H). The 
Coulomb energy ( |3.2|) is minimized by the state where (g(x)) = 0. It is realized when the 
electron density is homogeneous, 

(p(x)) = (<ft(x)tp(x)) = e 2 <^» (<pt{x)ip(x)) = po- (4.12) 
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In the semiclassical approximation we obtain 



N N 



H(pW) = He 2 ^>(^(x N ) • • ■^{x 2 )ip\x 1 )ip{x 1 )^{x 2 ) ■ ■■ip{x N )\$). (4.13) 

r=l r=l 

We insert a complete set \ n ){ n \ = 1 between two operators (p'(xi) and ip(xi). When the 
state |$) contains N electrons, only the vacuum term |0)(0| survives in the complete set 
because (p(x r ) decreases the electron number by one. Hence, iV-body wave function (|4.8|) is 
given by 

N 

%lx]=Po /2 U e ~ {AiXr)) > ( 414 ) 

r=l 

up to an irrelevant phase factor. On the other hand, to suppress the kinetic energy we 
impose the LLL condition ( |3.10|) , as requires $ v [x] to be analytic. Consequently, it follows 



that (A(x)) = constant or B c g = 0, which is possible only at v — 1/m from (|3.12|) . Namely, 
the ground state is realized only at the magic filling factor v — 1/m. At v — 1/m the 
ground-state wave function is given by $<p[x] = constant in terms of composite bosons, and 
therefore by the Laughlin wave function ( |2.2j ) in terms of electrons. In this way the Laughlin 
state is proved to be the ground state in the improved CB theory. 

For the sake of completeness, we briefly recall the results of the corresponding analysis 
with use of bare composite bosons. We can derive the equations similar to ( |4.10| ) and 
( [4.11] ), which verifies the mapping ( |3.7| ). The semiclassical ground state is similarly given 
by ^tj,(x) = constant in terms of bare composite bosons. The wave function in terms of 
electrons is singular, 

The state does not belong to the lowest Landau level. We can remove this singular short- 
distance behavior by including a higher order perturbation correction |2,15|. Nevertheless, 
it makes the naive CB theory less attractive. 



V. SEMICLASSICAL ANALYSIS 



We analyze excitations on the QH state. A priori two types of excitations are possible, 
that is, perturbative and nonperturbative ones in terms of the density fluctuation g(x) and 
its conjugate phase x( x )- To carry out a perturbative analysis, we parametrize the bare field 
as 4>{x) = e lx ^\/po + g(x) in terms of the density deviation g(x) and its canonical phase 
x(x)- The dressed field <p(x) is a nonlocal operator due to the factor e~ A ^ with ( |3.14| ), 

tp{x) = e-^gW*) Vpo + g( x ). (5.1) 

Substituting (|5.1|) into the Hamiltonian (|4.1|) and expanding various quantities in term of 
q{x) and x{ x )i we obtain the perturbative Hamiltonian, which is found to be identical to 
the one analyzed already in the bare CB theory ||, as should be the case. The result 
is that there exist no perturbative fluctuations confined to the lowest Landau level. We 
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conclude that all excitations in the lowest Landau level are nonperturbative objects. The 
improved theory confirms this assertion by showing explicitly how they are created on the 
ground state. Indeed, any excited state is represented as in ( |4.10| ), which is a nonlocal object 
because the creation operator (p*(x) is a nonlocal operator as in ( |5.1|) . 



A. Vortices 

We first examine excited states when the wave function ( |4.8| ) is factorizable, Q v [x] = 
uj\z] = Yl r u ( z r)- m this case we can easily make a semiclassical analysis of the one-point 
function (<f(x)) = u(z) by setting 



'■' ~ '"^^ ) Vpo+~Q^) =u(z), (5.2) 



-■A(x) ix(x) , 



based on the parametrization (53). Here and hereafter, we use the same symbols A(x), 
g(x) and x{ x ) a l so f° r the classical fields. When an analytic function u)(z) is given, (|5.2|) is 
an integral equation determining the density deviation g(x) confined to the lowest Landau 
level. 

We transform (|5.2[ ) into a differential equation. The Cauchy-Riemann equation for the 
analytic function ( |5.2|) yields, 

dj(A(x) - In a/po + q(x)) = -e jk d k x(x). (5.3) 
Using (|3.13| ) we obtain that 



-^V 2 ln(l + ^) _ e ( x ) = uQ v (x), (5.4) 
which we call the soliton equation, where 

Q V i x ) = ^-£jkdjd k x{x) = -^E jk djd k \Yiuj(z) (5.5) 

Z7T Z7T1 

is the topological charge density associated with the excitation. This is nonvanishing since 
lnuj(z) is a multivalued function unless uj(z) = constant. The topological current is 

Qu ( x ) = 7r £ nvxdvd x x(x), (5.6) 

with Q v (x) = Qq(x). It is a conserved quantity, d^QY h {x) = 0. Topological solitons are 
generated around the zeros of lu(z) according to the soliton equation (|5.4j) . The density 
modulation is induced in order to confine the excitation within the lowest Landau level. 
The topological charge is evaluated as 

Q V = J d 2 xQ v (x) = dx k d k \nu(z), (5.7) 

where the loop integration § is made to encircle the excitation at infinity (|a;| — > oo) provided 
oj{z) is regular everywhere. The topological charge is uniquely determined by the asymptotic 
behavior of the classical field <p(x), 
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<p(x) — > y/poZ q , as \z\ — > oo, (5.8) 
for which the electron number is 

AN = d 2 xg{x) = -vQv = -vq, (5.9) 



as follows from the soliton equation ( |5.4j) . The electric charge carried by this soliton is 
—eAN = uqe. It is a hole made in the condensate of composite bosons. We may as well 
derive the electric number (|5.9|) more directly from the parametrization (|5.2|). We take the 



asymptotic behavior \z\ — > oo in (|5.2| ) and equate it with 

g(x) — > 0, x( x ) ~ > Q&7 A(x) — > — 5 In \z\, as \z\ — > oo. (5.10) 

Taking the limit \x\ — > oo in ( [3.14J) we recover the quantization of the electron number ( p.9| ). 
Actually we have determined the coefficient y/po in the asymptotic behavior (|5.S| ) in this 
way. 

A topological excitation carries a quantized topological charge, which has to be created 
all at once. It cannot be created by an accumulation of perturbative fluctuations, as agrees 
with the perturbative result mentioned at the begining of this section. In the absence of the 
Coulomb term all these excitations are degenerate with the ground state, which explains 
the degeneracy in the lowest Landau level at v < 1. The degeneracy is removed since any 
density modulation acquires a Coulomb energy, 

(He) = 7T / d xd v-\ r = i v 1 ~o~i 5 - n 

2e J \x- y\ e£ B 

where 7 is a constant of order one. 

An explicit example is given by a vortex (quasihole) sitting at 1 = 0, whose wave 
function is $ v [x] = Y\ r z r up to a normalization factor, or uj{z) = ^/p^z. In this example 
the topological charge (|5.5| ) is concentrated at the vortex center, Qq(x) = S(x). A crude 
approximation of the soliton equation ( |5.4|) reads 

Y^ 2 Q( X ) - Q{x) = u6{x), (5.12) 

since v = 27cp £ 2 3 . Its exact solution is g(x) = —{u/-k£ 2 b )Kq{s) with s = \/2r/£B and 
Ko(s) the modified Bessel function. This is a rather poor approximation because of its 
singular behavior, g(x) — ► —00, at the vortex center (x — 0). A better approximation is 
given by 

2 

g(x) = -p (l + s-^-)e- s , (5.13) 



6 

which has the correct behavior both at x = and |x| 3> £b- Furthermore, it has the correct 
topological charge. For a numerical analysis, it is convenient to set p(x) = poe u ^ in the 
soliton equation (|5.4j), as yields JTB|, 

d 2 u 1 du , „ f ,i / _ 1 a\ 

— + -— + ! = e< s \ (5.14) 
as z s as 

for s > 0. The result of a numerical analysis shows that the density modulation is well 
approximated by (|5.13 ), as in Fig.|l|. The Coulomb energy is given by ([5.11 ) with 7 ~ 0.39. 
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VI. QUANTUM HALL FERROMAGNET 



We proceed to analyze the QH system with the SU(2) symmetry. The electron field 
i/j a (x) has the index a =f, L It denotes the electron spin in the monolayer QH system 
with the spin SU(2) symmetry, or the layer index in a certain bilayer QH system with the 
pseudospin SU(2) symmetry. For definiteness we analyze the monolayer spin system in this 
paper. 



A. Hamiltonian 

The Hamiltonian depends on the electron spin through the Zeeman energy term, 

H z = -9*^bB J d 2 xS z (x), (6.1) 

with S z = i( , 0'rt'0T_'0lt'0l) j where g* is the gyromagnetic factor and \xb the Bohr magneton. 
Each Landau level contains two energy levels with the one-particle gap energy g*nsB. The 
lowest Landau level is filled at v = 2. We consider the case where the Zeeman energy is 
much smaller than the Coulomb energy. Though one Landau level contains two degenerate 
energy levels in the vanishing limit of the Zeeman coupling (g* = 0), the system becomes 
incompressible at v = 1/m. The physical reason is the Coulomb exchange energy, as we 
now see. 

We define the bare CB field <j) a (x) and the dressed CB field ip a (x) by 

(f) a ( x ) = e- i9 ( x) i/> a (x), if a {x) = e- A{x) <j) a {x), (6.2) 

where the phase field Q(x) and the auxiliary field A(x) are given by ( |3.5| ) and ( |3.14| ), re- 
spectively, with the total electron density p(x) = J2 a ^p a ^ (x)^p a (x) = J2 a f aji (x)e 2A ^ip a (x). 
The Hamiltonian is 

H = ^mY.I ^r\x){P 2 x + P 2 y W{x) + H C + H Z (6.3) 
= ^^(*)) V^^(x) + ^ c + H C + H z , (6.4) 



with the Coulomb term He and the Zeeman term Hz- The Coulomb term depends on the 
deviation g(x) of the total electron density from the average density, g(x) = p(x) — p , as in 



B. Wave Function 

We may decompose the bare CB field into the U(l) field 4>(x) and the SU(2) field n a (x), 
<j) a {x) = <f){x)n a {x), ^n a \x)n a {x) = 1. (6.5) 
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The field n a (x) is the complex-projective (CP) field |13 |, whose overall phase has been 



removed and given to the U(l) field 4>(x). The spin operator is expressed as 

S a (x) = ip(x)S a (x), (6.6) 

where 

£*(*) = n\x)r a n{x\ n(x) = ( ^ ) . (6.7) 

In terms of the dressed CB field the decomposition reads 

^ a {x) = <p(x)n a (x), <p(x) = e- A{x) <f){x). (6.8) 



The SU(2) component n a (x) is common between the bare and dressed fields ( |6.5| ) and ( |6 . 8[ ) : 
It is a local field. On the other hand, (f(x) is a nonlocal field due to the factor e - - 4 ^- 1 as in 
the spinless theory. 

The ground state minimizes both the Coulomb and Zeeman energies. The Coulomb 
energy is minimized by the homogeneous electron density, (p{x)) = po- The Zeeman energy 
is minimized when all electrons are polarized into the positive z axis, (n*(x)) = 1 and 
(n^(x)) = 0. The ground state is unique, which we denote by \go). 

The two-component CB field is <&(x) = ip(x)n(x). With the Hamiltonian ( |6.4|) , the LLL 
condition for the state |$) is 

A$ (x) |$ ) = . (6.9) 

Because of this condition the iV-body wave function is analytic, 

%[x] = (0mx 1 )$(x 2 )---$(x N )\$) = Q[z], (6.10) 

where Q[z] is totally symmetric in iV variables. When it is factorizable, Q[z] = Y[ r ^{ z r), 
the electron wave function reads 



$ n=ii(^))^n- (6.11) 



The ground-state wave function is 



*N=n(J) $LN N> (6.i2) 



where all spins are polarized. 



C. Spin Texture 



We consider a spin texture given by performing an SU(2) transformation on the ground 
state \g ), 
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\$)=e i0 \g o ), (6.13) 

where is its generator, 

= E / dV»S a (x) = W d 2 Qf- q S a q . (6.14) 

a a 

The spin texture is described by the classical sigma field s a (x) = ($|S a (x)|$), which we 
parametrize as 



s x (x) = <t(x), s y (x) = \J\ - a 2 {x) sini?(x), s z (x) = a/1 - a 2 (x) cosf9(x). (6.15) 
It is classified by the Pontryagin number fl3|], Q p = J (1 2 xQq(x), with the topological current 

Qu( x ) = -^-^abce^\s a d u s b d x s c . (6.16) 

It is absolutely conserved, c^Q^ = 0. 

The spin texture ( |6.13| ) does not belong to the lowest Landau level. The excitation 
energy is naively given by (^/\H\^>) with the Hamiltonian ( |6.4j) , to which the kinetic term 
and the Zeeman term contribute but the Coulomb term does not. The kinetic energy is of 
the order of u c , as implies the spin stiffness of this order. Such an excitation is impossible 
at sufficiently low temperature. It is necessary to excite merely the component |$) of |$) 
belonging to the lowest Landau level. Furthermore, it is also necessary to extract the LLL 
component of the potential term since it kicks out the state out of the lowest Landau level. 

D. LLL Projection 

The LLL component |$) is extracted from the spin texture |$) by extracting the LLL 
component {{f a {x))) from the "wave packet" f a (x) in the generator of the SU(2) trans- 
formation ( |6.14| ). This turns out to replace the plane wave e lxq with [i7,IE| 



((e ixq )) = e-h H2 B& x i ( 6 . 17 ) 

in the Fourier representation of f a (x), where X = (X,Y) is the guiding center. We call 
{{e 1XQ )) the LLL projection of e lxq . The generator ( |6.14j ) is projected as 

O = E / ^((f a ^)))S a (x) = Y,j dVV*' ( 6 - 18 ) 



where 



Similarly we define 



S a = (2ti)- 1 / d 2 x((e- iqx ))S a (x). (6.19) 



p q = &*T X J d 2 x((e- iqx ))p(x) (6.20) 
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for the electron density operator. 

From the commutation relation [X, Y] = —i£% between the X and Y components of the 
guiding center, we obtain 

[e^,e ipX ] = 2ie i(9+p)x sinf^— ], (6.21) 

with qAp = q x p y — q y p x . The translation e iqx is Abelian, but the magnetic translation e iqX 
is non-Abelian. It governs the symmetric structure of the two-dimensional space after the 
LLL projection. It is straightforward to derive the following WooXSU(2) algebra [I8|, 



[Pp> Pq] = sin^l 2 ?] exp[-f pg], (6.22a) 

[S a P , p q \ = ^ p+q sin[£|f ] exp[f pg], (6.22b) 
fe^] = ^ abc S c P+q cos[£l^ ]exp[fpg] 

+ i-5 ab p p+q sin[£|f ] exp[f pq], (6.22c) 



based on the algebra ( 6.21|) of the magnetic translation. 



We make the LLL projection of the Hamiltonian (6.4), whose result we denote by H. 
The Coulomb term reads 

Hc = \J d 2 xd 2 yg(x)((V(x - y))) e {y) = vr j d 2 qV(q)p^ q p q , (6.23) 
where V(q) is the Fourier transformation of the potential V(x) = e 2 /e\x\. 



E. Exchange Energy 



We evaluate the energy by making a perturbative expansion of the spin texture 

around the ground state, 

H cS = ($\H |$) = (g \H\g ) - (g \[O, H)\g ) + ■■■. (6.24) 

Making a straightforward algebraic calculation, making a gradient expansion and taking the 
lowest order term in /* we obtain the exchange energy |T5|,P , P!9| from the Coulomb term 

(P3), 



(6.25) 



where s a (x) is the nonlinear sigma field with p s = ue 2 / (16^/2ne£B)- It describes the spin 
stiffness. Combining the exchange energy and the Zeeman energy we obtain the effective 
Hamiltonian, 



H cS = ^p s J2 J d 2 x[d k s a \ 



x)] 2 -^g*p B B j d 2 xs z (x). 



(6.26) 
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The exchange energy arises because of the following reason: The local spin rotation has 
components in higher Landau level since it is not a symmetry of the Hamiltonian. Only its 
LLL component is excited at sufficiently low temperature. Since the LLL components of 
the spin operators and the density operator do not commute as in (|6.22b|) , the local spin 
rotation induces a local density modulation and affects the Coulomb energy. 

The spin-stiffness term is precisely the nonlinear-sigma- model Hamiltonian ||. Though 
the term is derived perturbatively in the present framework, there is much numerical evi- 
dence |3, 20 , |21| that it captures correctly the long-distance physics associated with Skyrmion 
excitations. 



F. Goldstone Mode 

The Zeeman effect is quite small in actual samples. We consider the vanishing limit of 
the Zeeman term (g* = 0). According to the effective Hamiltonian ( |6.26D , the energy is 



minimized for any constant value of the sigma field, s(x) = So =constant. Hence, there 
exists a degeneracy in the ground states as indexed by s . The choice of a ground state 
implies a spontaneous magnetization, or a quantum Hall ferromagnetism. When a continuous 
symmetry is spontaneously broken, there should arise a gapless mode known as the Goldstone 
mode and quantum coherence develops spontaneously. 

When we include a small Zeeman effect, the ground state \go) is chosen where so = 
(0, 0, 1). We can treat the Zeeman interaction as a perturbation because it is much less im- 
portant than the Coulomb interaction. The key property of the QH ferromagnet is that it is 
a macroscopic coherent state though the coherent length is finite, where all spin components 
are simultaneously measurable, s a (x) = 2pQ 1 (S a (x)), with extremely good accuracy. An 
evidence is the existence of coherent excitations such as Skyrmions. 

The Goldstone mode describes small fluctuations of the CP field around the ground state 
( |6.12| ). Up to the lowest order of the perturbation in the CP field n(x), it is parametrized 



as 



n\x) = 1, nHx) = ^=i, (6.27) 



with [((x),(Ky)} = $( x ~v)- The LLL condition ( |6.9| ) yields two conditions, 

^(x)|$> = 0, A C(x) |$ ) = , (6.28) 

up to this order. Although they look similar, they describe very different excitation modes. 
As in the spinless QH system, ip(x) is a nonlocal field and generates extended objects. On 
the other hand, ((x) is a local field, and it describes the Goldstone mode. 
We may relate ((x) to the classical sigma field ( |6.15|) , 



((( x )} = y-^{a(x) + itf(x)}. (6.29) 



The effective Hamiltonian (|6.26| ) is recognized as a classical counterpart of the quantum 
version, 



15 



tf eff = ^ f d 2 x[d k (\x)d k ((x) + r L 2 C\x)C(x)], (6.30) 
Po J 



on the coherent state. Here, £^ is the coherent length, 



. 2ft (2 J r) 1/4 ^ . . 



where 



5 = 2/ . 6-32 

is the ratio of the Zeeman energy to the Coulomb energy. We call it the normalized g-factor. 
We have ~ in typical samples at B ~ 10 Tesla, where g ~ 0.02. In the momentum 
space the effective Hamiltonian reads 

H eS (k) = E k (l( k , (6.33) 
with [Cfc, Q] = $(k — I), and the dispersion relation is 

E k = —A; 2 + g*n B B. (6.34) 
Po 

The Goldstone mode has acquired a gap E = g*fisB. 



VII. TOPOLOGICAL EXCITATIONS 

We analyze topological excitations on the QH ferromagnet. We use the semiclassical 
approximation. When the iV-body wave function is factorizable, the one-point function is 
analytic, (ip a (x)) = u a (z). From ( |6.8|) , the one-point function is parametrized as 

e-^V^Vpo + Q{x)n a {x) = u a (z), (7.1) 

since |0(x)| 2 = Po + q(x). Here and hereafter, all fields are classical fields. When the wave 
function u a (z) is given, the electron density g(x) and the spin field S a (x) are determined 
by this equation. There are two types of excitations associated with the U(l) part and the 
SU(2) part of the CB field. The U(l) excitation has a characteristic length £b, while the 
SU(2) excitation has no scale provided the Zeeman term is neglected. 



A. Vortex Excitations 



The U(l) excitation is generated on the spin-polarized ground state (|6.12|) when d k x{x) ^ 



and d k n a (x) = in (|7.1|) . We may set (cp^) = 0. The one-point function (<p^(x)) is 
essentially Abelian, and the Cauchy-Riemann equation for (|7.1|) yields precisely the same 
soliton equation ( |5.4| ). The topological charge density is given by ( |5.5| ) with an analytic 
function u(z) = lj'(z). The U(l) excitation is the vortex. The Coulomb energy of the vortex 
excitation is given by ( |5.11| ) with 7 ~ 0.39. There is no antivortex excitation. Instead of it 
an electron is placed into the spin-down state, as would increase the Coulomb energy of the 
same order as the vortex excitation and the Zeeman energy. 
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B. Skyrmion Excitations 



The SU(2) excitation is generated on the spin-polarized ground state ( |6.12| ) when 



d k x(x) = and d k n a (x) ^ in (O). The CP field is 



n <*( x )= l £M (7 2) 



yielding the wave function 



^S kyrmi onN = n(^^)) ^NN- (7.3) 



A simplest choice is given by 

^SkyrmionN = ( j $LnN, (7.4) 

with a positive integer q, which describes a classical Skyrmion with scale n sitting at the 
origin of the system. It is clear in ( |7.4|) that the Skyrmion is reduced to the vortex in the 
limit k — > 0, where there is no distinction between the U(l) and SU(2) excitations. Because 
the vortex is regarded as the small limit of the skyrmion, we do not make a clear distinction 
between them in the QH ferromagnet. 

For the Skyrmion (|7.4| ) the classical sigma field ( |6.15| ) is calculated as 

s* = v/R^cos^), s y = -vT^fsm(q9), s z = g ~ { * B J* , (7.5) 

and the Pontryagin number density ( |2.11| ) as 

The spin flips at the Skyrmion center, s = (0,0,-1) at r = 0, while the spin-polarized 
ground state is approached away from it, s = (0, 0, 1) for r ^> k£ b . 

A Skyrmion excitation modulates not only the SU(2) part but also the U(l) part via the 
relation ( pM|) . The Cauchy-Riemann equation reads 

dj(A(x) - In a/po + Q{x)) = -e jk K k , (7.7) 

where 

K k = -iJ2n a *d k n a . (7.8) 



From (|7.7| ) the same soliton equation as ( [5.4D is derived 

^V 2 lnfl + ^ 

4vr V p 



-^V 2 ln( 1 + ^l) - g(x) = vQ p (x), (7.9) 
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but the topological charge density now reads 



Q p (x) = ^ £]k d 3 K k {x). (7.10) 
Ztx 



It is a straightforward calculation to show that the charge ( 7.10|) is identical to the 
Pontryagin number density ( |2.1l ). 



The topological charge is evaluated as 



Q = J d 2 xQ p (x) = ^-f dxjKj(x), (7.11) 

where the loop integration § is made to encircle the excitation at infinity (\x\ — > oo). The 
electron number associated with the topological soliton is 

AN = J d 2 xg(x) = -vQ p . (7.12) 

The topological charge is determined by the asymptotic value of the CP field (|7.2p. We 
find Q p = q for the Skyrmion ( [7.51) . The electron number of this soliton is AN = — vq: It 
represents the number of electrons removed by the Skyrmion excitation. 
The Skyrmion spin is given by 

AN S = ~J d 2 x{2S z (x) - Po } = 1 J d 2 x{p - p(x)s z (x)}. (7.13) 



It seems to diverge logarithmically for the Skyrmion ( |7. 5| ) with q = 1. This is a fake since 
the the Zeeman term breaks the spin SU(2) symmetry explicitly and introduces a coherent 
length £l into the SU(2) component. The Skyrmion configuration ( |7.5| ) is valid only within 
the coherent domain because the coherent behavior of the spin texture is lost outside it. By 
cutting the upper limit of the integration at r = k£l/2 in ( [7.13j ), we obtain 



where £x is the coherent length given by ( |6.31| ) and g is the normalized g-factor ( |6.32| ). 

The density modulation around the Skyrmion is governed by the soliton equation ( [7.9|) . 
This equation has formally the same expression as the soliton equation (|5.4j) for the vortex 



excitation. We may obtain an approximate solution in the two limits, the large Skyrmion 
limit (k, ^> 1) and the small Skyrmion limit (k < 1). First, in the large limit we may solve 
([7.91) iteratively as 



v 2 



q(x) = -vQ p {x) - —V 2 Q p (x) + ■■■. (7.15) 
87rpo 



We may approximate it as 



g(x) ~ -vQ p {x) = - { l_ B f for « » 1, (7.16) 

7r [r 2 + {tu^YY 



18 



for the Skyrmion with q = 1, where we have used (|7.6| ). It agrees with the formula due 
to Sondhi et al. ||. On the other hand, the topological charge Q p (x) is localized in the 
small limit, Q p {x) — > qS(x) as k — > in ( [7.61 ). Hence, the solution is given by the vortex 
configuration, 

£(*)^- Po (l + ^-^)e-^, for «:«!, (7.17) 



which has been derived in (|5.13|) . 



We evaluate the excitation energy of a Skyrmion with q = 1 at v = 1. It consists of the 
exchange energy, the electrostatic term and the Zeeman term, where the exchange energy 
( |6.25| ) is exactly calculable for one Skyrmion, 

^Skyrmion = ^ P s + ^ / d^j/ f^^ + 9 *fl B BAN a , (7.18) 

2e J \x-y\ 



where the Skyrmion spin AN S is given by ( [7. 13| ) . It is calculated by using ( [7.1 6| ) and ( |7. 14| ) 
for a large Skyrmion, 



,2 



E, 



Skyrmion 



(3 gn ( V2n 



si 



B 




+ — In + 1 

32 k 2 V 32g 



(7.19) 



with (3 = 37r 2 /64. For a sufficiently small Skyrmion, the Coulomb energy is calculated with 
the vortex configuration and is given by ( |5.11| ). 

The Coulomb energy increases for a smaller Skyrmion while the Zeeman energy increases 
for a larger Skyrmion. The optimized scale k is obtained by minimizing the total energy 

(El). 



The Skyrmion excitation energy is given by ([7.19 ) together with ( 7.20 ). 



The parameter (3 describes a strength of the Coulomb energy. In general it is a function 
of the size k, obeying f3(n) —>■ 0.39k as k — > and /3(k) — ► 37r 2 /64 as k — > oo. We expand it 
around an arbitrary point as /3(k) = /3q + /3i(k — k ). The /t-dependence of the Skyrmion 
energy is the same as in ( |7.19| ) with /3 = /3o — /3iKq. The optimized scale is given by ( |7.20| ) 
with the same replacement. The parameter (3 is calculable only if the soliton equation ( [7.9| ) 
is solved for a Skyrmion wih arbitrary scale k. Furthermore, there will be a correction 
to it from a finite thickness of the layer. In this paper we treat (3 as a phenomenological 
parameter. 

The creation energy of a Skyrmion-antiSkyrmion pair will be given by 2i?sk y rmion if they 
are sufficiently far apart one another. However, it is not this quantity that is observed 
experimentally. The activation energy A is usually determined by the Arrhenius formula, 

R xx oc ex P[-77^L ( 7 -21) 

where R xx is the magnetoresistivity. Various factors affect the activation energy [0. For 
instance, impurities make Coulomb potentials around them and lead to a Landau-level 
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broadening. Phenomenologically, their effects result in a subtraction of a certain amount of 
offset r offset from the Skyrmion creation energy, 



A = 2E. 



Skyrmion 



offset ) 



(7.22) 



where r offset increases with disorder. In FIG.@ we have fitted the experimental data due to 
Schmeller et al. |]TT| by the formula ( |7.19p , where we have used j3 = 0.24 and an appropriate 
offset r offset phenomenologically for each curve. The theoretical curve reproduces all the data 
remarkably well. The Skyrmion spin is estimated that AN S ~ 3.7 at ~g — 0.015 (B = 3.05 
Tesla). This estimation is consistent with the Hartree-Fock result || and the experimental 
result based on Knight-shift measurements IJlO 



VIII. DISCUSSIONS 

We have studied QH ferromagnets based on the improved CB theory. We have investi- 
gated excitations confined to the lowest Landau level. Charged excitations are Skyrmions, 
carrying the quantized charge ue at the filling factor v — 1/m. Small Skyrmions are identi- 
fied with vortices. We have successfully explained the experimental data due to Schmeller 



et al. 11 
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FIG. 1. The density moduration around a vort ex with q = 1 is plotted. The solid curve is 
obtained by solving the differen tial equation ( p. 14 ) numerically. The dashed curve is drawn by 
using the approximate formula ( |5.13| ). 
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FIG. 2. Theoretical curves versus experimental data due to Schmeller et al. [11] for the acti- 
vation energ y in QH ferromagnets. All the data are fitted excellently by the theoretical formula 
( 7.22j ) with ( [T.19D , where an appropriate offset r o fj set is assumed for each sample. 
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